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A theoretical model for understanding the formation of interfacial thin films is presented, 
which combines density functional theory calculations for interface energies with thermody- 
namic modeling techniques for multicomponent bulk systems. The theory is applied to thin 
film formation in VC-doped WC-Co cemented carbides. It is predicted that ultrathin VC films 
may exist in WC/Co interfaces at the high temperature sintering conditions where most of 
the WC grain growth occurs, which provides an explanation of the grain growth inhibiting 
effect of VC additions in the WC-Co system. 
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1. Introduction 

Ultrathin films are essential for many applications in nanoscale science and technol- 
ogy. Considerable scientific interest is devoted to thin oxide films at metal surfaces 
[IH3], for which applications are found in catalysis, semiconductor devices as well as 
mechanical wear and corrosion protection. Correctly assessing the thermodynamic 
stability of such films is an important issue New thin film surface alloys have 
been discovered and the stability of nano-sized patterns in the ultrathin films in- 
vestigated with respect to strain and chemical effects [H [6] . Formation of ultrathin 
films is also observed in metal-ceramic materials after sintering of nanoparticles 
into a bulk nanocrystalline material [7J. 

In this Letter, we consider an interface between two dissimilar materials, a metal- 
ceramic interface, and the propensity to form, at the interface, an ultrathin film 
of a secondary phase, whose corresponding bulk phase is thermodynamically un- 
stable. We develop a model for the thin film formation and stability as function of 
temperature and show that it can be stabilized by interfacial effects. 

For modeling of interfaces with complex bonding such as metal-ceramic inter- 
faces, first-principles calculations in the framework of density functional theory 
(DFT) has proved to be a very useful tool [8HT0]. However, phase stabilities at 
high temperatures are cumbersome to predict with sufficient accuracy using DFT 
and therefore we combine DFT calculations with driving forces of nucleation 
from well- developed thermodynamic modeling techniques [12]. These modeling 
techniques describe bulk phases but lack information on interfacial effects. 
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We apply the developed methodology to WC-Co cemented carbides. These ma- 
terials are produced by means of powder metallurgy, where powders of carbides 
and metal Co are sintered together into a hard and dense material with excellent 
mechanical properties. The materials are used in cutting and wear resistant tools 
and are of considerable industrial importance. Much effort has been invested in de- 
veloping nanostructured cemented carbides by reducing the mean WC grain size in 
the finished product, because a fine and homogeneous microstructure has superior 
mechanical properties [13]. To control the microstructure and retain a fine grain 
size throughout the sintering, dopants in the form of VC, Cr3C2, TiC, TaC or NbC 
powders are added, of which VC has proved to be the most effective inhibitor in 
ultrafine alloys [Hj. The grain growth inhibiting effect of the dopants is well estab- 
lished, whereas the inhibiting mechanism on the atomic scale is not understood. 
Experimental high-resolution electron microscopy (HREM) studies [15H19] have 
shown the existence of V- and Cr-rich films of cubic carbide structure at WC/Co 
interfaces in VC- and Cr3C2-doped materials, respectively, at low temperatures, 
after liquid phase sintering and cooling. In order to predict the possible growth 
inhibiting effect of such films, it is decisive to determine whether these can exist at 
high temperature liquid phase sintering conditions, where most of the grain growth 
occurs and where HREM cannot be applied. 



2. Model 

2.1. General description 

The interface energy 7 can be expressed as the excess Gibbs energy per unit inter- 
face area due to the presence of an interface [20| , 




G is the Gibbs energy for the entire considered system, which contains an inter- 
face of area A. is the number of atoms in the system of component i with a 
corresponding chemical potential ^ determined by the reservoirs with which the 
system is in equilibrium. In the present case, the reservoirs are stoichiometric WC 
and a binder phase. The latter consists of mainly Co (and will be denoted "Co") 
but also of dissolved W, C, and M atoms, where M = V, Cr, Ti, etc. depending on 
the doping. 

We consider the propensity to form a thin film consisting of N layers of carbide 
phase, here denoted MC, at the interface between WC and the binder phase. For 
sufficiently thick films, the interface can be regarded as split into a WC/MC in- 
terface and a MC/Co interface separated by a slab of MC. The energy for such a 
film-covered interface can then be decomposed as 

Tfiim = 7wc/mc + 7mc/Co + N (A#mc + eMc) , (2) 

where 7wc/mc (7mc/Co) is the interface energy between WC and MC (MC and 
Co) . The last term of Eq. ([2]) represents the energetic cost per layer of building the 
MC phase and, thus, increases linearly with the number of layers N of the film. For 
the corresponding energy per formula unit we use the notation A(7mc = SAgMC, 
where S is the interface area of one formula unit. The energy A^mc is the negative 
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Figure 1. A schematic illustration of 7fli m as function of film thickness N. The circles denote the energy 
in an atomistic treatment and the solid line follows the decomposition in Eq. |2|. 

driving force of nucleation. Assuming a stoichiometric carbide phase, it equals 

= 9MC ~ MM - fJ>c, (3) 

where gyic is the Gibbs energy per formula unit of MC and /%[ (nc) is the chemical 
potential of M (C). A positive value of A$mc thus implies that the MC phase is 
unstable. We consider coherently strained thin films and the last quantity of Eq. ([2]), 
CMC; is the strain energy per layer needed to bring the film into epitaxy with the 
underlying WC. 

The energy 7gi m should be compared with the energy 7 W c /Co °f a WC / Co inter- 
face not covered by any MC film. The film is stable if 7 mm < 7wc/Co an d, hence, 
a positive value of 

A7MC = 7wc/Co - (7WC/MC + 7MC/Co) ( 4 ) 

indicates that film formation is possible. An upper limit of the film thickness -/V"ij m j t 
is found where the energy of the film-covered interface matches the energy of the 
uncovered interface, 7 nmi = 7wc/Coi which yields 

AT A7MC / R v 

AWt = x= — • ( 5 ) 

&9mc + CMC 

If iViimit is larger than a few atomic layers, the interface may contain a thin film. 

The decomposition of 7 nmi in Eq. ^ is valid only if the film is sufficiently thick, 
so that the interior of the film is bulk- like and the interface energies are well-defined. 
For thin films, 7mm will deviate from the linear dependence on N, which implies 
that 7fii m can have a minimum as function of TV" for a finite film thickness. In equi- 
librium, this optimal film thickness N eq will be attained. A schematic illustration 
of 7fiim is given in Figure [T] 



2.2. Specific model and method 

In the present paper we concentrate on VC-doped WC-Co. Results for other MC 
dopants will be presented elsewhere [21] . WC has a simple hexagonal lattice struc- 
ture. After liquid-phase sintering in Co, the WC grains are predominantly delimited 
by their basal (0001) and prismatic (1010) surfaces |22| . Here, we focus on the basal 
surface, which can be either W- or C-terminated. We choose the WC surface and 
the cubic VC phase to be oriented according to 



(0001) wc II (lll) vc , [1210] wc || [110] vc . (6) 
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This orientation relationship between WC and MC has been experimentally estab- 
lished for both V- and Cr-doped WC-Co [MES]. 

Interface calculations are performed using DFT and the current computational 
methodology has been developed in a series of papers [24r:26]. Computational re- 
sults have been compared with experimental data |24| [25] and the use of the gen- 
eralized gradient approximation (GGA) for the exchange-correlation effects within 
DFT in interface calculations has been motivated [26]. Further computational de- 
tails can be found in [27J. 

To model the interfaces, we use standard supercell slab geometries, where the 
interfaces are made coherent by straining the VC and Co phases to match the 
WC basal surface. The associated strain energies are minimized by relaxing the 
interplanar distances to account for the Poisson effect. In the calculation of interface 
energies, the strain energies are canceled, so that 7wc/VC) 7vc/Co> an d 7wc/Co 
represent the chemical contribution to the interface energy and do not depend on 
system size. 

To define an absolute value of the interface energy, a definite value of the carbon 
chemical potential fic is required. Within the two-phase region (WC+Co) of the 
W-C-Co system, [Iq can be varied in a small interval limited by the formation of 
either graphite or a metal-carbide phase (W,Co) 6 C [2S]. We assume carbon to be 
in equilibrium with graphite and use the corresponding computed value for (jlq [lOj. 



3. Results and discussion 
3.1. Interface energies 

In the (0001) wc and (lll) vc directions, the WC and VC phases consist of al- 
ternating layers of metal and carbon atoms. Assigning the positions A to (0,0), 
B to (3,!) and C to (§,3) with respect to hexagonal lattice vectors in the 
(0001) wc plane, the stacking sequence for WC is . . . W^CbWaCb • • • and for VC 
. . . VaCcVbCaVcCb • • • along [0001] wc and [lll] vc , respectively. We have calcu- 
lated interface energies for all possible stacking sequences involving the A, B and C 
positions and varying atomic terminations of the phases. Only stoichiometric VC 
compounds are considered. To assign a value of 7wc/VC) we Pi°k the lowest energy 
among the different stackings and terminations. The resulting interface energy is 
found to be slightly negative, 7wc/vc = —0.03 J/m 2 j2"7] . 

The predicted WC/VC interface structure can be compared with experimen- 
tal findings. By comparing experimental and simulated HREM images of a 
WC(0001)/VC(111) interface in a VC-doped WC-Co, a particular interfacial stack- 
ing has been suggested [29]. We find that the lowest energy is obtained for 
. . . aCbV aQcV bC>a ■ ■ ■ [2Z] which is in complete agreement with the stacking 
suggested in Ref. [29J . 

The WC/Co and VC/Co interfaces need to be treated with more concern, since 
at relevant sintering temperatures, the Co binder is in its liquid phase. A solid- 
liquid interface energy is usually considered as stemming from two separate parts 
(see, e.g. [30]): a chemical contribution due to differing atomic bonding at the 
interface, and a structural contribution due to disturbance of the liquid's structure 
caused by the interface. One should notice, that only the difference between the 
two interface energies 7wc /Co an d 7vc /Co enters into the expression Q for the key 
quantity A7vc- Therefore, the geometrical similarity of the WC(0001) and VC(lll) 
surfaces implies that structural contributions to A7VC will to large extent cancel. 

We choose to model the WC/Co and VC/Co interfaces in the same geom- 
etry as the WC/VC interface, which means that we are effectively modeling 
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Figure 2. The optimal atomic structure of a WC/VC/Co interface containing two layers of 
stoichiometric VC. Top: Atomic positions projected onto a plane with normal [1210] wc- Bottom: The 
positions of the V, C and Co atoms closest to the interface projected onto the interface plane. The six 
depicted atomic configurations of Co are used to calculate a mean value of the carbide/Co interface 
energy for the thin film configuration. The same set of relative translations and rotations of Co with 
respect to the topmost carbide layer is also used in calculating WC/Co and VC/Co interface energies. 

WC(0001)/Co(lll) and VC(lll)/Co(lll) interfaces with different Co stackings. 
Following the preceding argument, we do not expect the particular choice of Co 
orientation to influence A7VC in any decisive way, although existing experimental 
observations of this orientation between WC and solid Co [31, 32] do motivate 
its use as model geometry. At high temperatures, the motion of liquid Co atoms 
will constantly change the atomic structure at the interface. Therefore, we use a 
mean value of energies for different stackings of Co in the interface to approximate 
7wc/Co an d 7vc/Co) respectively. In Figure [2] we depict the six different atomic ar- 
rangements of Co used in calculating the mean energy of the carbide/Co interface. 
The resulting value for the interface energy can be regarded as a good approxi- 
mation for an incoherent interface [33]. The termination of the WC or VC slabs 
(metallic or C) is consistently chosen to yield the lowest mean interface energy. 

We find that for WC/Co, the interface is C-terminated with an energy 7wc/Co 
1 .01 J/m 2 [27]. For VC/Co, the interface is V-terminated with 7vc/Co = 0-00 J/m 2 
|27j . The resulting difference in interface energy is Atvc = 1-03 J/m 2 , which we 
henceforth will assume to be temperature independent. 

The final result of A7VC does not strongly depend on the choice of averaging 
of the various Co stackings in the interfaces. For instance, by following the above 
procedure but disregarding the two most high-energetic configurations in each in- 
terface, A7VC increases by 0.14 J/m 2 . The impact on A7VC due to ordering of 
the liquid Co atoms in the vicinity of the WC or VC interfaces is therefore ex- 
pected to be maximally of this size. Temperature-dependent deviations from our 
estimation of A7VC may arise from differences in the vibrational spectrum of V, 
W, and C atoms at the interface as compared with their respective reference states. 
Following the procedure of Reuter and Scheffler [34] , we estimate this effect to be 
maximally ±0.3 J/m 2 . Considering both these effects, we expect A7VC to be deter- 
mined within an accuracy of ±0.35 J/m 2 . By accounting for a mixed carbide layer 
instead of the stoichiometric VC compounds considered so far, further entropic 
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3.2. Thin film energy and composition 

In addition to A7VC; the energetic costs, Agvc and evc> to build the VC thin films 
are needed to predict stability. To compress VC from its calculated equilibrium 
lattice parameter in the interface, ayc/v^ = 2.944 A, to the calculated WC lattice 
parameter, awe = 2.919 A, an energy of only eye = 0.011 J/m 2 (0.005 eV per VC 
formula unit) is required. 

The driving force of nucleation, — Agvc> depends sensitively on the temperature 
and composition of the system. The most stable cubic carbide is not stoichiometric 
VC, but a mixed (V,W,Co)C x phase. In its at 1593 K experimentally measured 
composition, V constitute 75% and W 23% of the metallic atoms and the C vacancy 
concentration is 11% [28J. To determine the driving force of nucleation we have used 
the Thermocalc software [12] with a database for Cr and V additions in the W- 
C-Co system |28j. The database has been assessed from experimental observations 
using the Calphad method and models temperature dependent Gibbs energies of 
the phases occurring in the W-C-Co system. 

In Figure [3l the negative driving force of nucleation Ag( V ,w,Co)c for a 
(V,W,Co)C 3 , phase of optimal composition with respect to Gibbs energy in a W- 
C-Co-V system in equilibrium with WC and graphite is given as a function of tem- 
perature for varying V/Co atomic ratio y». For low temperatures, Ao( V ,w,Co)c = 
for all considered doping levels, meaning that stable (V,W,Co)C x precipitates exist 
at these chemical conditions. The sharp increase of A<?(v,w,Co)C a t temperatures 
around 1550 K is due to the melting of the binder phase, which significantly in- 
creases its solubility of V. At the doping levels yi and 1/3 of Figure [3j the thermody- 
namic calculation (which lacks any model for interface stabilization) predicts that 
all (V,W,Co)C a . phase would dissolve as the binder melts. A discontinuous change 
in the slope of Ao( V ,w,Co)c is a l so found near 1300 K for y\. This change is as- 
sociated with the transition of the solid Co binder from its ferro- to paramagnetic 
phase. We conclude that for relevant doping concentrations, Ao( V ,w,Co)c varies 
between and 0.1 eV (corresponding to and 0.22 J/m 2 per layer of (V,W,Co)C a; ). 

We can now summarize our results for the propensity of thin film formation. 
In Figure [5j the solid line corresponds to the result from Eq. ([2]) for VC. Here, 
we set Agvc = 0.06 eV, which is a reasonable choice based on the results from 
our thermodynamic modeling at relevant sintering temperatures. With this value 
of AovCj Eq. ([5]) suggests that stable films of thickness less than iViimit = 7 m ay 
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exist. 

In the interface calculations, we have approximated the film carbide phase with 
a stoichiometric VC phase. The real film may well have a composition similar to 
the bulk (V,W,Co)C a . phase, though finding the optimal composition would require 
numerous large interface calculations. To assess the effect of composition, we have 
performed a calculation for a mixed (V0.75, Wo.2s)C phase, whose composition is 
close to the experimentally measured one of bulk (V,W,Co)C a . [28]. We obtain 
7(Vo.75,Wo.25)C/Co = -0.15 J/m 2 < 7vc/Co [27] and, therefore, we expect that for 
the true composition, the interfacial effect will further stabilize the thin film. 



3.3. Explicit atomistic layer- by -layer treatment 

The expression for 7mm in Eq. ([2]) is only valid if the film is sufficiently thick. For 
thinner films, the details of the atomic layer-by-layer structure has to be taken into 
account. We have therefore explicitly modeled the atomic layer- by-layer structure 
of thin films of cubic VC at the WC/Co interface. The number of possible stacking 
sequences of VC in the WC/VC/Co interface is large, but can be decreased by 
omitting obviously high-energetic configurations. In total, we have considered 76 
different VC stacking sequences with six corresponding Co stackings for each of 
these [27] . 

In Figure |4j the energy of a WC / Co interface containing a VC film is given as 
a function of the number of metallic V layers in the film for different VC stacking 
sequences. The results for the stackings corresponding to the minimal energy in 
the WC/VC interface described earlier are circled. In the plot, we show the result 
without adding any contribution from the energetic cost per layer to build the 
thin film (Agvc + eve = 0)- The result should then converge to 7wc/vc + 7vc/Co 
for sufficiently thick films. We find that the result is close to the asymptotic value 
already for two V layers, while the result for one V layer deviates both energetically 
and structurally with . . . CbW^CbVcCo . . . being the optimal stacking in this case. 
The optimal structure in the case of two V layers is depicted in Figure [2j 

Our final results are presented as circles in Figure [5] where we again set 
A<?vc = 0.06 eV. To calculate 7fiim(A0; we take for each number of V layers 
iV the minimum energy among the stackings and terminations presented in Fig- 
ure [4] and add N (Agvc + eve)- The fast convergence of the interface energies 
toward the asymptotic value of Eq. ([2]) implies that a thin film with two V lay- 
ers will be the most stable configuration for a rather wide interval of Agyc 
0.01 eV < A<7vc ^ 0.2 eV. For no value of Agvc will one V layer be the most 
stable configuration. With the values used in Figure [5] we find that the interface 
with a VC film is favored by 0.68 J/m 2 compared to the interface without segre- 
gated VC. This value lies well outside our estimated error bound of ±0.35 J/m 2 
indicating that in equilibrium, the WC(0001) surface is indeed covered by VC. 

Our prediction of a few stable V layers agrees very well with HREM imaging 
[7J [181 EE [22] , where thin films on WC(0001) in VC-doped material are seen to have 
a thickness of two and sometimes a few additional metallic layers. An interpretation 
of the additional layers is that they precipitate as Agvc — > during cooling. The 
VC-covered surface is also a likely nucleation site for further growth of the VC 
phase, which sometimes builds up as a precipitate on WC(0001) during cooling 

P21E9]. 
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Figure 4. Interface energy 7fli m as a function of the number of V layers TV in a cubic VC film in the 
WC/Co interface for different stackings of VC. + (x) denotes a C-terminated (metal-terminated) 
carbide/Co interface. The stacking sequence that corresponds to the minimal WC/VC interface energy is 
circled. No energetic cost per layer to build the thin film is added (Agvc + g vc = 0). 




Figure 5. Interface energy 7fli m as a function of the number of V layers TV in a cubic VC film in the 
WC/Co interface. Agyc is se t to 0.06 eV per formula unit, which is an appropriate value for relevant 
doping conditions and sintering temperatures. A thin film of two V layers will be the most stable 
configuration for a rather wide interval of AgvCi 0.01 eV < Agvc $ 0.2 eV. 

3.4. Grain growth inhibiting effect 

It is generally believed that the rate limiting step of grain growth in WC-Co systems 
is the dissolution and/or reprecipitation reactions occuring at WC/Co interfaces. 
Several authors have therefore attributed the growth inhibiting effect of V to a 
proposed interference with these reactions either through blocking active growth 
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centers at the WC surface [351 ES] or through the presence of segregated V at 
the WC surface [HI [19] , which would hinder W diffusion into or out of the WC 
grain. Our prediction that ultrathin VC films may be present and stable also at 
liquid-phase sintering temperatures is supportive of the latter explanation. 



4. Conclusions 

We have developed an atomistic model for understanding thin film formation at 
interfaces in bulk composite materials. We combine extensive density functional 
theory (DFT) calculations for interfaces and thin films with thermodynamic mod- 
eling techniques for multicomponent systems. The latter technique is used to obtain 
relevant chemical potentials for the various components as function of temperature 
and doping conditions. 

The theory is applied to thin film formation in VC-doped WC-Co cemented 
carbides. It is predicted that under realistic doping conditions, ultrathin films may 
exist at high temperature liquid phase sintering conditions where most of the grain 
growth occurs. 

By explicit DFT modeling of the thin films as function of the number of atomic 
layers, we find that the most stable configuration consists of only two V layers over 
a large temperature range and we show how the atomistic picture approaches the 
thick film limit with increasing number of layers. 

The films, whose corresponding bulk phase is thermodynamically unstable, are 
stabilized by interfacial effects. Compared with the original interface in the absence 
of a thin film, there is an energy gain associated with the creation of the two sub- 
interfaces involved in the thin film formation, which more than compensates for 
the energetic cost of the thin film itself. 

In conclusion, our results explain the formation of VC films in WC/Co interfaces. 
The presented method can be applied to theoretically assess the stability of other 
carbide films in WC/Co interfaces. 
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COMPUTATIONAL DETAILS 

We employ DFT as implemented in the Vienna ab- 
initio simulation package (VASP) [lj. The exchange- 
correlation functional is approximated in the Perdew- 
Burke-Ernzerhof (PBE) scheme [2]. The plane- wave 
pseudopotential method with projector augmented wave 
(PAW) [3, 4] potentials is used. The plane-wave energy 
cutoff is set to 400 eV in all calculations. Partial oc- 
cupancies are set with the method of Methfessel-Paxton 
[5] of first order with a smearing parameter of 0.05 eV. 
Atomic relaxations are performed until all atomic forces 
are smaller than 0.02 eV/A. 

We consistently use non-spinpolarized calculations, 
which effectively means that all considered phases are 
treated as paramagnetic. Although ferromagnetic at low 
temperatures, Co is paramagnetic above its Curie tem- 
perature of 1396 K [6]. The Curie temperature is low- 
ered by solutes, and for Co in equilibrium with WC and 
graphite, it is around 1324 K (from the model of [6]), 
which is below the temperatures of concern in the cur- 
rent paper. 



INTERFACE CALCULATIONS 

In the interface calculations, the periodic supercell con- 
tains two interfacing slabs of the corresponding phases 
and at least 14 A of vacuum. This setup allows calcula- 
tion of energies for individual stacking sequences in the 
WC (0001) /VC (111) interface, which is not possible in a 
WC/VC slab setup without vacuum, since the two inter- 
faces formed in such model would be inequivalent due to 
the ABC stacking of the VC phase. The carbide phases 
are modeled as slabs of 5 + 4 atomic layers and the Co 
phase as a slab of 8 atomic layers. Analogous geometries 
are used for interfaces containing films of a secondary 
phase. 

For the WC and VC compounds, the smallest possible 
cell size (lxl atoms in the (111) or (0001) planes) is used, 
where the Brillouin zone is sampled with a T-centered 
grid with a division of 13 x 13 x 1 points. For the cal- 
culation of 7(v 75 .w 25)c/Co> a ce U of 2 x 2 atoms in the 
(111) plane is used with a fc-point sampling of 7 x 7 x 1. 
Each (111) plane contains 3 V and 1 W atoms, and the 
placement of W atoms in successive (111) planes is chosen 
so that W-W interatomic distances are maximized. For 
bulk calculations, fc-point sampling is done consistently 
with the interface calculations. Considering all technical 



errors, interface energies are hereby converged to a max- 
imal error of 0.015 eV per interfacial unit corresponding 
to 0.03 J/m 2 . 

The energies described in Tables I, II, III and IV are 
taken as total energy differences from DFT calculations. 
As an example, 

7wc/Co =-j (E - N W fi w - A C Mc - Nco/ico) ^ 
- ctco - owe, 

where E is the total energy of the interface system con- 
taining Nyj (Ac, Aco) W (C, Co) atoms, /i\y, Mc and 
/iCo are calculated in separate bulk calculations. Due 
to our choice of geometrical model, which contains vac- 
uum, the surface energies of WC and Co, <r\yc and ctq , 
obtained in separate surface calculations are also sub- 
tracted. All calculations are done in the same strain 
state corresponding to fixed lattice parameters in the 
(0001) wc || (Hl) vc || (Hl)co P lane at the calculated 
equilibrium lattice parameter of WC. Thus, the strain 
energy in VC and Co is canceled in the calculation of 
interface energies. 

When constructing a WC/Co interface containing a 
VC film, several stacking sequences are possible for each 
number of layers of VC. From the values in Table I, we 
conclude that all stacking sequences that contain two 
subsequent layers of C or two layers of atoms placed on 
top of each other (e. g. ... | Va • • •) can be disre- 
garded due to their high energy. Furthermore, we only 
study stacking sequences that are compatible with a nacl 
stacking. The remaining number of stacking sequences is 
still rather large. Counting the number of possible stack- 
ing sequences for one metallic layer in the film, we find 
16 different configurations. For two and more metallic 
layers, 20 VC stackings are considered. For every VC 
film stacking, six different stackings of Co on top of the 
WC+VC slab are modeled in order to give a mean value 
for the metal-ceramic interface energy. Up to four metal- 
lic layers in the film have been explicitly modeled, which 
gives a total of 456 DFT calculations. 
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TABLE I: Interface energy for different stackings in the 
WC (0001) /VC (111) interface. Energies are given in J/m 2 . 



Stackin; 


I sequence 


7 


...WaCbWa 


VaCcVbCa • • ■ 


3.79 


...WaCbWa 


VaCbVcCa • • ■ 


3.87 


...WaCbWa 


VbCcVaCa • • ■ 


1.10 


...WaCbWa 


VbCaVcCb... 


1.27 


...WaCbWa 


VcCbVaCc • • • 


1.32 


...WaCbWa 


VcCaVbCc • • • 


1.49 


...WaCbWa 


CaVcCbVa • • . 


5.20 


...WaCbWa 


CaVbCcVa • • ■ 


5.19 


...WaCbWa 


CcVbCaVc... 


1.38 


...WaCbWa 


CcVaCbVc • • • 


0.92 


...CbWaCb 


VaCcVbCa • ■ ■ 


-0.03 


...CbWaCb 


VaCbVcCa • ■ ■ 


0.39 


■ . . CbWaCb 


VbCcVaCb • • • 


3.04 


. . . CbWaCb 


VbCaVcCb • • • 


3.03 


...CbWaCb 


VcCbVaCc • • • 


0.83 


...CbWaCb 


VcCaVbCc • • • 


0.29 


...CbWaCb 


CaVcCbVa • ■ ■ 


7.81 


...CbWaCb 


CaVbCcVa . . . 


6.28 


. . . CbWaCb 


CbVcCaVb • • • 


4.13 


. . . CbWaCb 


CbVaCcVb • • • 


4.09 


...CbWaCb 


CcVbCaVc • • • 


8.22 


...CbWaCb 


CcVaCbVc • • • 


9.72 



TABLE II: Interface energy for different stackings in the 
WC (0001) /Co (111) interface. Energies are given in J/m 2 . 



Stacking sequence 



7 



...WaCbWa 


CoaCobCocCoa . . • 


2.59 


...WaCbWa 


CoaCocCobCoa . . • 


2.53 


...WaCbWa 


CobCoaCocCob • • • 


0.07 


...WaCbWa 


CobCocCoaCob • • • 


0.35 


...WaCbWa 


CocCoaCobCoc • • • 


0.64 


...WaCbWa 


CocCobCoaCoc • • • 


0.94 




Mean 


1.19 


...CbWaCb 


CoaCobCocCoa ■ • • 


-0.02 


...CbWaCb 


CoaCocCobCoa ■ • • 


0.22 


. . . CbWaCb 


CobCoaCocCob • • • 


3.27 


. . . CbWaCb 


CobCocCoaCob • • • 


3.16 


...CbWaCb 


CocCoaCobCoc • • • 


-0.21 


. • • CbWaCb 


CocCobCoaCoc . . . 


-0.36 




Mean 


1.01 



TABLE III: Interface energy for different stackings in the 
VC (111) /Co (111) interface. Energies are given in J/m 2 . 



Stackin 
•CaVbCcVa | 
•CaVbCcVa | 
•CaVbCcVa | 
•CaVbCcVa | 
• CaVbCcVa | 
• CaVbCcVa | 



g sequence 
CoaCobCocCoa • 
CoaCocCobCoa • 
CobCoaCocCob • 
CobCocCoaCob • 
CocCoaCobCoc ■ 
CocCobCoaCoc ■ 





Mean 


0.00 


. . . VbCcVaCb 


| CoaCobCocCoa • • • 


-0.99 


. . . VbCcVaCb 


| CoaCocCobCoa • • • 


-0.73 


. . . VbCcVaCb 


| CobCoaCocCob ■ • • 


3.45 


. ■ . VbCcVaCb 


| CobCocCoaCob ■ • • 


3.27 


. . . VbCcVaCb 


| CocCoaCobCoc . . . 


-1.53 


. . . VbCcVaCb 


| CocCobCoaCoc . . . 


-1.71 




Mean 


0.29 



7 
0.62 
0.68 
-0.51 
-0.24 
-0.36 
-0.16 



TABLE IV: Interface energy for different stackings in the 
(V0.75, Wo.25)C (111) /Co (111) interface. Energies are given 
in J/m 2 . 

Stacking sequence 7 



...Ca(VW)b C c(V,W) ;1 1 


CoaCobCocCoa • • • 


0.64 


...CU(V,W) B Cc(V,W) A 1 


CoaCocCobCoa • • • 


0.71 


...Ca(V,W) b Cc(V,W)a 1 


CobCoaCocCob . . • 


-0.72 


...Ca(V,W) b Cc(V,W)a 1 


CobCocCoaCob . . . 


-0.45 


...Ca(V,W) b Cc(V,W) a 1 


CocCoaCobCoc • . . 


-0.63 


...Ca(V,W) b Cc(V,W) a 1 


CocCobCoaCoc . . . 


-0.42 


Mean 




-0.15 


...(V,W) B Cc(V,W) A C B 1 


CoaCobCocCoa ■ ■ • 


-1.03 


...(V,W) B Cc(V,W) A C B 1 


CoaCocCobCoa . . . 


-0.78 


...(V,W) b Cc(V,W) a Cb 1 


CobCoaCocCob • • • 


3.25 


...(V,W) b Cc(V,W) a Cb 1 


CobCocCoaCob • • • 


3.09 


...(V,W) b Cc(V,W)aCb 1 


CocCoaCobCoc . . . 


-1.55 


...(V,W) b Cc(V,W)aCb 1 


CocCobCoaCoc . . . 


-1.72 


Mean 




0.21 



